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ABSTRACT 


Coherent  light  beams  propagating  through  the  atmosphere  undergo 
considerable  phase  perturbations  due  to  fluctuating  temperature  structures  in 
the  atmosphere.  Understanding  and  measuring  these  structures  on  a  real-time 
altitude  dependent  basis  is  inherent  to  successful  deployment  of  ground  based 
lasers  and  particle  beams. 

One  method  used  to  detect  these  temperature  structures  is  an  acoustic 
profiler,  or  echosounder.  Of  immediate  interest  is  the  ability  of  high  frequency  (5 
kHz)  planar  array  echosounders  to  rapidly  detect  low  level  turbulence  (below 
200  meters)  and  quantify  the  results. 

This  thesis  involves  design  improvements  of  previously  developed 
echosounder  arrays  and  associated  software.  Particulary,  this  thesis 
demonstrates  that  tighter  packing  of  elements  in  a  planar  acoustic  array 
produces  better  side  lobe  reduction  than  less  densely  packed  arrays.  This 
results  in  higher  energy  density  in  the  main  lobe  and  increased  performance. 

Also  included  in  this  thesis  is  a  computer  method  which  allows  relatively 
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Optical  degradation  of  coherent  light  in  the  atmosphere  is  strongly  a 
function  of  small  scale  temperature  gradients  which  change  the  local  index  of 
refraction  in  the  atmosphere  [Ref.  1].  These  localized  structures  in  the 
atmosphere  fluctuate  rapidly,  changing  in  size,  density,  location,  velocity  and 
interaction  with  neighboring  structures.  There  is  also  considerable  variability 
with  diurnal  temperature  changes,  annual  weather  variations  and  regional 
climatological  differences.  Understanding  the  mechanisms  of  this  turbulence 
and  measuring  it  on  a  real-time  local  level  are  intrinsic  to  high  power  earth 
based  lasers. 

These  small  scale  turbulence  structures  in  the  atmosphere  play  an  even 
greater  role  with  respect  to  acoustic  absorption  and  reflection  [Ref.  2].  It  is  well 
known  that  acoustic  pulses,  on  the  order  of  1  to  6  kHz  frequency,  are  ideal  in 
defining  turbulence  structures  in  the  first  few  kilometers  of  the  atmosphere,  with 
lower  frequencies  having  greater  ranges.  Until  recently,  the  most  popular  type 
of  acoustic  sounder  was  a  monopole  speaker/receiver  mounted  in  a  parabolic 
dish  [Refs.  3  and  4],  Several  problems  with  this  method  include  difficulty  of 
transportation,  high  amplitude  side  lobes  at  large  angles  and  poor  altitude 
resolution  due  to  the  low  frequencies  involved.  These  deficiencies  can  be 
overcome  by  means  of  planar  arrays  [Refs.  5  and  6],  and  this  thesis  was 
designed  to  amplify  this  point. 
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Working  under  Walters,  Weingartner  and  Wroblewski  [Refs.  7  and  8]  again 
demonstrated  the  feasibility  of  this  technique  by  building  and  testing  a  25 
element  5  kHz  planar  array.  The  result  was  an  extremely  effective  low  altitude 
echosounder  which  was  less  expensive,  more  portable  and  had  better  side 
lobe  reduction  than  previous  monople  sounders. 

During  the  course  of  this  previous  work  several  questions  arose 
concerning  the  planar  array.  If  the  elements  were  more  densely  packed,  then 
could  the  size  of  the  side  lobes  be  significantly  reduced,  thus  increasing  the 
accuracy  of  the  beam?  Could  the  array  be  modeled  on  a  computer  to  produce 
the  best  design?  Could  the  enclosure  be  improved? 

This  thesis  is  an  attempt  to  answer  these  questions,  and  involves  the 
development  of  a  dense  packed  5  kHz  hexagonal  array  and  associated  sound 
deadened  enclosure.  Additionally,  various  computer  programs  using  fast 
Fourier  transform  algorithms  were  studied  as  a  first  approximation  in 
diagnosing  acoustic  side  lobe  suppression. 

The  results  indicate  that  densely  packed  planar  acoustic  echosounders 
produce  better  side  lobe  suppression  than  more  open  arrays  and  are  accurate 
tools  in  the  identification  of  low  level  turbulence. 
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The  atmosphere  induces  random  phase  perturbations  on  an 
optical  beam.  Diffraction  of  these  small  scale  phase  variations  produces  a 
larger  beam,  restricting  the  ability  to  focus  energy  on  a  target.  Optical  phase 
perturbations  in  the  atmosphere  stem  from  temperature  fluctuations  carried  by 
the  turbulent  velocity  field.  These  temperature  fluctuations  vary  in  size, 
cascading  from  tens  of  meters  in  diameter  down  to  centimeters  and  smaller. 
These  structures  appear  to  be  most  variable  and  inhomogeneous  at  the 
interface  between  stratified  layers,  particularly  near  the  earth’s  surface. 

Within  the  first  few  hundred  meters  of  the  surface,  these  temperature 
structures  are  dominated  by  the  surface  heat  flux,  resulting  from  differing  rates 
of  heat  transfer  between  the  earth  and  the  atmosphere.  During  the  day  the 
surface  heats  much  more  rapidly  and  non-uniformly  than  the  atmosphere.  An 
insulating  blanket  of  air  near  the  surface  heats  and  rises  off  the  ground,  more  or 
less  in  sheets.  Warmer  sections  rise  more  rapidly,  quickly  disrupting  the 
uniformity  of  these  sheets.  Cooler  air  moves  in  to  replace  the  rising  air,  further 
disrupting  the  sheets  and  causing  vertical  thermal  plumes  to  form.  These 
thermal  plumes  continue  upwards  in  the  atmosphere,  slowly  mixing  and 
increasing  the  temperature  of  the  lower  atmosphere.  This  cycle  is  continuous, 
with  the  actual  transfer  of  energy  occurring  by  small  scale  molecular  diffusion. 
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At  night  the  cycle  reverses,  with  the  surface  cooling  off  much  more  rapidly 
than  the  air  via  Planck  radiation.  In  this  case,  as  the  air  gives  up  energy  to  the 
cooler  surface,  it  produces  a  turbulent  boundary  layer  (via  velocity  interaction 
with  a  stationary  surface)  which  remains  relatively  close  to  the  ground.  This 
cooling  process  is  generally  less  turbulent  than  the  heating  process,  due  to  the 
reduced  heat  flux  at  the  boundary,  although  intense  turbulence  continues  to 
exist  at  the  interface  of  the  horizontally  stratified  layers. 

There  is  also  a  neutral  event,  during  the  transition  between  night  and  day, 
where  the  surface  temperature  equals  the  ambient  air  temperature.  During  this 
time  the  net  heat  flux  between  the  air  and  the  surface  changes  sign  and  passes 
through  a  minimum.  The  atmosphere  generally  becomes  optically  stable  during 
these  brief  periods. 

The  structural  variability  is  clearly  evident  in  optical  instruments  designed 
to  measure  the  turbulence  of  the  atmosphere  [Refs.  10  and  11],  but  altitude 
correlation  with  optical  instruments  alone  is  difficult.  Acoustic  echosounders 
alleviate  this  difficulty. 

Acoustic  echosounders  are  devices  which  generate  and  propagate 
energy  into  the  atmosphere  and  measure  the  backscatter  return  from  turbulent 
structures  as  a  function  of  altitude.  The  power  returned  from  these  turbulent 
structures,  Pr,  is  proportional  to  the  power  transmitted,  Pj.  The  equation 
relating  these  terms  has  been  summarized  by  Neff  [Ref.  12],  based  upon  work 
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by  Tatarski  [Ref  1],  and  is  repeated  in  equation  (1)  below. 

(1)  PR  =  PT  Er  Et  (ct/2)  (exp(-23R)}  (AG/R2)  p[R,/]  where, 

-  PR  is  the  electrical  power  returned  from  range  R, 

-  PT  is  the  electrical  power  transmitted  at  frequency  /, 

-  Er  is  the  efficiency  of  conversion  from  acoustic  to  electric  power, 

-  Et  is  the  efficiency  of  conversion  from  electric  to  acoustic  power, 

-  c  is  the  local  speed  of  sound  (m/s), 

- 1  is  the  acoustic  pulse  length  (s), 

-  A  is  the  antenna  area  (m2)  , 

-  R  is  the  range  to  the  scattering  volume  (m) 

-  exp  (-2I3R)  is  the  round  trip  power  loss  due  to  attenuation  where  3  is  the 
average  attenuation  to  and  from  the  scattering  volume, 

-  G  is  the  effective  aperture  factor  for  the  antenna, 

-  p  [r,/]  is  the  scattering  cross-section  per  unit  volume;  that  is,  the 
fraction  of  incident  power  backscattered  per  unit  distance  into  unit  solid 
angle  at  frequency  /. 


The  backscattering  cross  section,  p,  is  the  unknown  quantity  of  interest  in 
Eqn.  (1).  Tatarski  [Ref.  1]  uses  p  in  equation  (2)  (below)  to  find  CT2,  the 
temperature  structure  parameter. 

(2)  p[R./]  =  0.0039  k 1/3  (Cj2  /  T2)  where, 

-  k  is  the  acoustic  wave  number  =  2k  /  L, 

-  L  is  the  acoustic  wavelength, 

-  T  is  the  absolute  temperature  and 

-  Cj2  is  the  temperature  structure  parameter. 


The  temperature  structure  parameter,  Cj2,  is  the  mean  squared 


temperature  difference  between  separate  points  in  space, 

(3)  Ct2  =  {<T2(R2)-T1(Ri)>2}/{R2-R1)^  , 

and  provides  the  scattering  centers  which  produce  backscattered  energy  from 
the  turbulence.  <AT>2  empirically  increases  by  AR2/3.  Therefore,  it  is  normalized 
by  AR2/3  to  produce  a  relatively  constant  temperature  structure  parameter. 

Since  p  is  proportional  to  the  backscattered  return  power,  PR,  and  to 
the  temperature  structure  parameter,  CT2,  then  PR  is  proportional  to  C-p2. 
Solving  equation  (1)  for  p.  and  equating  equations  (1)  and  (2),  the  following 
expression  for  Cj2  is  given  for  an  acoustic  echosounder  of  wavelength  L; 

(4)  Cj2  =  1/0.0039  (1/ERET)T2k-1/3  (2/ct)(1 /AG)  PR/PT  {R2  exp(2RR)}. 

The  electric  to  acoustic  efficiency  factors  for  the  transducers  (speakers) 
were  measured  in  Reference  7  and  are  inherent  to  the  design  of  the  speakers. 
ER  =  Ej  =  0.5,  an  excellent  value  for  speaker  efficiency,  which  accounts  for 
their  use.  The  aperture  area,  A,  for  the  19  element  hexagonal  array  is  0.0866 
m2. 

The  effective  aperture  factor  for  the  array,  G,  accounts  for  the  non-uniform 
antenna  directivity.  G  is  a  difficult  quantity  to  calculate  or  measure.  It  represents 
a  normalized  average  of  the  acoustic  energy  within  the  transmit/receive 
antenna  beam  lobe  and  is  clearly  summarized  by  Probert-Jones  [Ref.  13]. 
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Based  upon  Reference  13  and  studies  by  Hall  and  Wescott  [Ref.  14],  a  value  of 
G  =  0.4  is  currently  used  in  computations.  This  quantity  is  presently  undergoing 
extensive  research  by  Dr.  Walters. 

Another  difficult  factor  in  Equation  (4)  involves  the  average  atmospheric 
attenuation  loss  factor,  6,  which  is  strpngly  dependent  upon  frequency  and  the 
absolute  water  vapor  in  the  atmosphere  [Refs.  2  and  15].  This  dependence  is 
illustrated  in  Figure  1,  which  is  taken  directly  from  Reference  15.  The  water 
vapor  pressure  is  measured  via  a  digital  humidity/temperature  indicator,  model 
5165-A  from  Weathermeasure  [Appendix  A]. 
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Fig.  1.  Attenuation  versus  Water-vapor  Pressure  (mb.) 
for  Frequencies  around  1.6  KHz. 


Combining  these  known  values  into  a  constant  Q,  equation  (4)  becomes; 

(5)  Cj2  =  (Q)  T2^173  {PR  /  (ct)(  PT)}  R2  exp(2BR)  where, 

Q  =  {2Pr}/{0.0039  Pj  Er  Ej  A  G}. 

Wroblewski  [Ref.  8]  utilizes  this  in  the  software  for  the  echosounder  to 
calculate  a  15  minute  time  averaged  value  for  the  temperature  structure 
parameter  as  a  function  of  altitude.  Uncertainties  in  the  antenna  beam  shape, 
side  lobes,  rejection,  electronic  gains,  transducer  efficiencies,  electronic 
antenna  gains  and  the  atmospheric  attenuation  produce  uncertainties  in  the 
absolute  calibration  of  an  echosounder. 

Independent  verification  of  the  CT2  values  measured  and  calculated  by 
the  echosounder  are  necessary  to  confirm  the  accuracy  of  this  technique. 
Without  this  confirmation  it  is  still  possible  to  guage  relative  variations  in  the  low 
level  turbulence,  simply  with  the  echosounder  profiles.  These  profiles  are 
presented  and  discussed  in  the  following  chapters. 


Much  of  the  equipment  used  in  this  project  was  borrowed  from 
the  25  element  square  array  designed  by  Weingartner  and 
Wroblewski.  The  equipment  basically  involves  a  Hewlett  Packard 
(HP)  217  computer  for  data  generation  and  analysis,  an  HP  331 4A 
function  generator  (controlled  by  the  software)  which  produces  the 
initial  pulse,  an  amplifier,  a  switch/pre-amp  which  passes  the 
initial  pulse  then  amplifies  and  transfers  the  return  signal  to  the 
computer,  the  array  and  shroud  which  act  as  the  physical 
transmitter/receiver  and  a  filter  used  to  eliminate  stray  noise. 
These  units  will  be  broken  down  and  individually  discussed  in  this 
section.  The  software  developed  for  the  echosounder  by  Wroblewski 
was  not  significantly  altered  and  is  not  included  in  this  document. 
[Refs.  7  and  8] 

A.  ELEMENTS 

The  elements  of  the  array  consist  of  19  piezoelectric 
Motorola  brand  speakers,  model  #  KSN  1005A,  with  attached  horns 
of  model  #  KSN  1032A.  Ceramic  piezoelectric  drivers  were  chosen 


because  of  their  small  size,  excellent  (see  Figure  2  )  frequency 
response  at  5  kHz,  and  a  shorter  ’ring'  time  with  respect  to 
electrodynamic  speakers.  Any  speaker  acts  as  a  damped 
harmonic  oscillator  with  a  'ring'  time  (recovery  time)  inversely 
proportional  to  the  damping  constant.  The  damping  is  much  stronger 
in  Piezoelectric  speakers  than  in  conventional  electrodynamic 


elements. 


Typical  Frequency  Response 
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Figure  2.  Speaker  Frequency  Response 


An  initial  set  of  40  speakers  was  tested  to  find  19  which 
were  similarly  matched  in  frequency  response.  This  testing 

was  done  by  comparing  the  output  of  a  chosen  reference 
speaker  with  the  acoustic-to-electrical  conversion  of  each 
test  speaker.  By  clamping  the  reference  speaker  face-to-face 
with  the  test  speaker,  a  known  voltage  to  the  reference  should 
produce  a  range  of  voltages  from  the  test  specimens.  These  voltages 
were  compared  using  Lissajous  plots  on  a  Nicolet  Oscilloscope.  It 
was  more  important  to  obtain  a  set  of  identical  speakers  rather  than 
the  best  speakers,  in  order  to  achieve  a  more  uniform  final  array. 

The  horns  on  each  driver,  originally  circular  with  a  square 
flange  face  were  trimmed  to  a  circular  diameter  of  7.62cm  such  that 
they  fitted  close  packed.  Ideally  a  diameter  of  d=  L/2 
(l=wavelength)  would  provide  optimal  phase  spacing  of  elements, 
but  this  required  either  smaller  speakers  (  d  »  3.4  cm  )  which 
were  unavailable,  or  a  lower  frequency  (=2230  Hz),  which  was 
below  the  useable  range  of  available  speakers. 

Another  source  of  speaker  was  identified  as  a  Philips  brand 
piezoelectric  2002  EX  [Ref-  16]  with  better  frequency  response  at 


3kHz  than  the  Motorola  speakers  in  use.  Unfortunately  these  Philips 
speakers  are  no  longer  manufactured. 

B.  ARRAY 

The  speaker  elements  were  configured  in  a  close-packed 
planar  array  on  an  aluminum-covered  balsa  wood  board  to  absorb 
vibrations  (Fig.  3).  The  elements  were  wired  in  parallel,  reducing 
the  array  impedance  to  300  Q  /  19  =  15.8  Q.  The  mounting  box  was 
locally  constructed  of  sheet  aluminum  and  enclosed  a  layer  of 
semi-porous  foam  loaded  with  a  thin  lead  sheet  to  reduce  stray 
noise  signals.  Input/Output  was  via  a  coaxial  connector.  Bracketts 
were  attached  to  allow  360°  rotational  testing  from  2  axes  in  the 
anechoic  chamber. 


Figure  3.  Array  Dimensions 


C.  ENCLOSURE 


The  enclosure  was  designed  to  be  easily  reproducible, 
inexpensive,  portable,  and  easily  assembled  [the  earlier  enclosure 
took  some  time  to  set  up  and  break  down].  Optimally,  a  long 
cylindrical  enclosure  was  desired  to  reduce  side  lobe  scatter  and 
return  without  obstructung  the  main  lobe.  The  enclosure  chosen  to 
optimize  the  above  requirements  was  a  Rubbermaid  55  gallon  plastic 
trash  container,  with  a  height  of  84  cm,  having  a  slight  upward 
taper;  the  barrel  has  a  lower  inside  diameter  of  53  cm,  opening  to  a 
diameter  of  60  cm.  This  provided  a  finished  open  half  angle  of  about 
20°,  8°  more  than  optimally  desired  according  to  tests  in  the 
anechoic  chamber  with  the  unencumbered  array  [see  RESULTS'!.  A 
second  possibility  was  a  highway  barrier  barrel  [Appendix  A],  but 
these  were  much  more  expensive  and  involved  similar  dimension 
problems. 

The  inside  surface  of  the  enclosure  was  covered  with  lead  lined 
semi  porous  foam,  with  7.62  cm  convoluted  (egg  crate  shaped)  foam 
on  one  side,  and  2  cm  flat  foam  on  the  backing  side.  The  entire  foam 
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liner  was  attached  with  standard  contact  cement.  The  lead  lining 
amounted  to  a  density  of  about  1  Ib/sqft  which  is  sufficient  to 
reduce  one  way  noise  penetration  by  40dB.  The  convoluted  foam 
helped  break  up  internal  and  off  axis  scatter  (spectral  reflection) 
and  aided  in  side  lobe  reduction.  The  foam  used  turned  out  to  be  the 

last  of  a  salvage  batch  of  foam,  and  a  second  source  of  foam  was 

desired.  An  ideal  substitute  for  lead  lined  foam  (which  was 
difficult  to  obtain)  was  found  to  be  a  barium  impregnated  vinyl  foam 
(a  thin  rubbery  vinyl  layer  loaded  with  barium  and  sandwiched 
between  semi-porous  foam)  with  a  similar  density  of  1  Ib/sqft  and 
virtually  identical  sound  absorption  capability  (40dB).  Two 
sources  of  this  foam  were  identified  [Appendix  A], 

D.  ELECTRONIC  HARDWARE 

The  data  processing  system  consisted  of  an  HP  Model  #217 
computer  with  a  20  megabyte  hard  drive  and  2  megabyte  memory,  an 
INFOTEK  BC  203  Basic  compiler,  analog  to  digital  (A  to  D)  converter, 

INFOTEK  FP  210  floating  point  accelerator  (to  enhance  the  speed  of 

execution),  HP  9133  floppy  disk  drive,  monitor  and  HP  ink-jet 
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printer  (see  Figure  4).  The  HP  217  was  the  workhorse  of  the  system, 
triggering  the  function  generator,  receiving  and  processing  the 
return  pulse,  and  outputting  the  results  in  graphical  mode.  Without 
this  unit,  timely  acquisition  and  reduction  of  data  would  be 
impossible. 

The  HP  3314A  function  generator,  adjusted  and  triggered 
automatically  by  the  HP  217  computer,  was  used  to  generate  a 
specific  transmitted  signal.  This  signal  was  generally  100  cycles  of 
a  sinusodial  waveform  of  1.5  volts  (rms)  amplitude. 

A  QSC  model  1700  audio  amplifier  was  used  to  increase  the 
voltage  to  the  array  from  1.5  volts  (out  of  HP  3314A)  to  30  volts 
(rms).  This  model  was  chosen  because  of  its  ability  to  drive  a 
highly  reactive  load  (the  array).  A  lighter,  more  compact  amplifier 
could  easily  replace  this  item  of  equipment,  and  will  be  considered 
during  future  upgrades. 

A  pre-amplifier  (pre-amp)  (Fig.  5)  was  used  to  isolate  the 
transmitted  pulse  from  the  received  pulse  and  to  amplify  the  return 
signal  by  a  gain  of  approximately  11,000.  The  pre-amplifer  was 
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Fig.  4  Echosounder  System  Set-up 
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designed  and  built  by  Walters  and  utilizes  two  low  noise  solid  state 


operational  amplifiers,  an  LT  1007  and  LT  1037.  All  four  isolation 
rectifiers  were  of  type  IN  4000. 

The  important  properties  to  consider  for  an  effective  pre-amp 
were  complete  isolation  of  the  circuit  during  transmission  (to 
prevent  overloading  components)  and  immediate  signal  retrieval 
following  the  'ring'  time  of  the  array.  The  more  rapid  the  recovery 

time  of  the  array  (a  function  of  driver  capacitance  (Ca)  and  the  loop 
resistance  (R-j )),  the  quicker  the  pre-amp  can  respond  to  returned 
signals  and  the  smaller  the  surface  blind  zone.  For  this  array,  the 
capacitance  of  each  driver  was  Cq  »  0.1  p  F.  Thus  the  time  constant 


for  the  array  was 

RC  =  (C[))(#drivers)(Ri )  =  time  constant 
-{O.ricr6  F)(1 9)(270  Q)  =  0.513  msec. 

For  a  given  array,  resistor  R-|  provides  the  time  constant  of  the 
circuit. 

After  5  time  constants  there  is  virtually  no  'ringing’  left  on  the 
array,  and  this  sets  the  minimum  height  of  the  blind  zone.  Thus: 
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5  x  T  =  2.57  ms  and  Hmjn  =  (340  m/s)(5.32*10'3s)  «  0.87  meters. 

Therefore  the  blind  zone  stands  at  approximately  1  meter  above  the 

array,  which  is  not  noticeable  on  the  echosounder  printout. 

A  Rockland  Wavetek  model  #753A  filter  (115  dB/octave}  was 

used  to  isolate  noise  and  further  amplify  the  return  signal  prior  to 

computer  analysis.  The  input  side  of  the  filter  was  set  up  as  a  low 

pass  filter  at  5  kHz  and  gain  setting  of  O  dB  (xl);  the  output  side 

was  set  up  as  a  high  pass  filter  at  5  kHz  with  a  gain  setting  of  20  dB 

(xl 0).  This  produced  a  band  pass  filter  centered  at  5  kHz  with  a 

band  width  of  about  100  Hz  at  the  -3  dB  level.  The  over  all  gain  of 

the  filter  in  this  mode  was  measured  as 

output  685  m  volt 

-  =  -  =  9.67  -  10 

input  70.88  m  volt 

The  combination  of  the  pre-amp  and  the  filter  produced  an  overall 
return  signal  gain  of  (pre-amp  gain)(filter  gain)  =  (1 1094)(9.67) 

-  107,000. 

The  gain  of  the  system  is  thus:  107,000  *  270  Q/(  270  +  15.8  Q  ) 


IV. 


A.  LOBE  PATTERNS 


The  array  was  tested  in  the  anechoic  chamber  utilizing  a 
program  designed  by  Butler  [Ref.  17].  This  program  utilizes  the 
set-up  shown  in  Figure  6  and  allows  the  measurement  of  beam 
patterns  for  various  arrays.  A  band  pass  filter  was  used  to  increase 
gain  and  reduce  60  Hz  noise  from  the  equipment.  Present  on  all  the 
lobe  patterns  is  a  spurious  lobe  at  =  135°.  This  was  determined  to 
be  a  reflection  from  the  inside  door  handle  of  the  anechoic  chamber 
and  should  be  ignored.  All  lobe  patterns  have  been  normalized  to  the 
maximum  voltage  detected,  allowing  easy  comparison  between 
figures. 

The  hexagonal  array  was  rotated  around  two  axes  to  obtain  the 
maximum  sidelobe  field  of  the  array  (Figs.  7  and  8).  The  pattern  for 
the  square  array  is  shown  for  only  one  axis  of  orientation  (Fig.  9). 
Comparison  between  the  two  arrays  (Figs.  7  and  9)  shows  a  definite 
decrease  in  side  lobe  number  and  intensity  for  the  hexagonal  array 
versus  the  square  array. 
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While  the  array  was  being  tested  in  the  anechoic  chamber  a 
demonstration  of  its  detecting  ability  was  desired.  Using  a  heat 
gun,  the  array  was  set  up  in  the  chamber  exactly  as  it  would  be  in 
the  field,  except  a  Nicolet  model  #3091  digital  Oscilloscope  was 
used  to  detect  the  return  signal. 

The  heat  gun  was  tested  at  three  different  settings;  first  was 
'heat  and  air',  then  'air  only'  and  last  ’heat  only'.  With  the  heat  and 
air  both  on,  a  definite  fluctuating  signal  could  be  detected  on  the 
oscilloscope  screen.  With  'air  only'  there  was  virtually  no 
detectable  signal  and  with  ’heat  only*  the  signal  was  almost  back  to 
its  level  of  ’heat  and  air'.  No  measurements  were  taken,  but  the 
demonstration  was  repeated  several  times.  This  directly  illustrated 
that  the  turbulence  structures  in  the  atmosphere  are  mainly  a  result 
of  temperature  fluctuations  and  not  due  to  turbulent  velocity 
fluctuations. 

The  next  step  involved  testing  the  side  lobe  suppression  of  the 
enshrouded  array.  j  mechanics  of  this  operation  were 

considerably  more  i  rolved  than  simply  suspending  the  array  from  a 
rotating  motor.  Being  too  heavy  to  suspend,  the  shroud  was  taped  to 


a  rotating  turntable  and  centered  beneath  the  rotator  motor.  The 
motor  was  connected  to  the  shroud  via  aluminum  rods  and  clamps. 


When  the  rotating  motor  turned,  so  did  the  shroud,  allowing  the 
enshrouded  beam  pattern  to  be  measured  and  recorded  with  the  same 
software  as  used  previously. 

The  initial  beam  pattern,  (Fig.  10)  showed  that  the  shroud,  while 
significantly  reducing  side  lobes,  was  not  long  enough  to  completely 
eliminate  them.  Some  of  the  side  lobes  were  getting  through. 

A  first-cut  step  in  further  reducing  these  side  lobes  involved 
placing  a  15  cm  wide  sleeve  of  barium  impregnated  smooth  foam 
around  the  upper  inside  circumference  of  the  shroud.  The  result  (Fig. 
11)  was  not  promising,  and  demonstrated  that  further  refinement 
was  required.  One  thing  Figure  11  demonstrated,  however,  was  that 
specular  reflection  from  a  smooth  surface  increases  side  lobe 
activity  and  thus  convoluted  foam  shrouds  (vice  smooth  foam 
shrouds)  are  required. 

The  next  step  involved  placing  a  15  cm  wide  annular  ring  coated 
with  sound  deadening  foam  over  the  exit  of  the  shroud  (similar  to  a 
trash  can  lid  with  the  center  removed).  This  left  an  open  circle  of 
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48.25  cm  diameter  for  the  acoustic  waves  to  propagate.  This 
adjustment  slightly  decreased  the  side  lobes  (==  3  dB),  but  not 
significantly  (Fig.  12).  Further  narrowing  of  the  neck  via  more 
annular  rings  only  succeded  in  reducing  the  main  lobe,  with 
little  subsequent  change  in  side  lobe  activity. 


B.  FIELD  MEASUREMENTS 


Field  measurements  were  conducted  on  several  occassions 


over  a  six  month  period  at  various  locations  around  California.  They 


generally  coincided  with  optical  measurements  being  conducted  by 


Walters  [Refs.  9  and  10]  and  were  used  to  supplement  the  optical 


data.  Of  greatest  interest  were  mountaintops  above  the  marine 


boundary  layer  (=4000  ft.),  where  the  turbulence  should  be  a 


minimum. 


Two  mountaintops  in  particular  were  visited,  one  being  Lick 


Observatory  on  Mt.  Hamilton  and  the  other  being  Mt.  Wilson 


Observatory  at  Mt.  Wilson.  These  were  chosen  because  of  the 


availability  of  telescope  domes,  which  were  required  to  reduce  wind 


vibration  on  the  optical  components  and  protect  the  computers  from 
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dust  and  debris.  The  echosounder  itself  is  impervious  to  wind,  and 
can  be  operated  wherever  electricity  is  available. 

Lick  Observatory  was  visited  on  1  and  2  August  1987  and  some 
of  the  results  are  presented  in  the  following  figures.  All  times  are 
given  as  universal  time.  Subtract  seven  hours  to  obtain  local 
(daylight  savings)  time  at  the  site.  Both  days  were  similar  with 
respect  to  turbulence.  Therefore,  only  profiles  from  2  August  will 
be  presented  to  illustrate  the  ability  of  the  echosounder  to  detect 
and  present  the  rapidly  varying  nature  of  the  atmosphere  accurately. 
The  CT2  profiles  are  currently  unverified  and  should  be  ignored. 

Figure  13  was  taken  early  in  the  evening  and  demonstrates  the 
presence  of  convective  pluming.  This  pluming  vanishes  gradually, 
eventually  yielding  to  the  nighttime  structure  which  starts  in  Figure 
14.  By  this  time  a  definite  layer  has  set  in  and  surface  pluming  is 
no  longer  dominant.  Later  in  the  evening  a  strong  turbulence  layer 
descendes  (Fig.  15)  which  quickly  settles  [Fig.  16]  and  is  quite  dense. 
This  layer  remains  close  to  the  ground  the  rest  of  the  evening,  with 
light  turbulence  detectable  above  it. 


RANGE  IN  METERS  RANGE  IN  METERS 


Early  in  the  morning  (Fig.  17)  another  strong  layer  is  noted.  As 
the  sun  begins  to  rise  (Fig.  18),  commencement  of  the  daily 


convective  pluming  becomes  clearly  evident.  These  figures 


illustrate  the  dramatic  ability  of  the  high  frequency  echosounder  to 


profile  low  level  turbulence  in  the  atmosphere.  Interestingly, 


whenever  a  turbulence  structure  was  detected  on  the  profiler,  there 


was  also  noticeable  degradation  to  both  the  optical  instruments  in 


The  next  occasion  to  use  the  echosounder  occurred  on  6  and  7 


November  1987  with  a  field  trip  to  Mt.  Wilson  Observatory. 


Unfortunatley,  no  echosounder  data  could  be  obtained  due  to  severe 


radio  frequency  [RF]  interference  from  approximately  10  radio/TV 


transmitters  co-located  at  the  observatory.  The  trip  was  not  a  loss, 


however,  as  it  was  used  to  correct  a  subtle  flaw  in  the  design  of  the 


pre-amp. 


The  problem  was  traced  to  the  'floating  point  ground'  design 


of  the  pre-amp.  The  coaxial  cables  leading  to  the  pre-amp  made 


ideal  antennas  at  the  200  megahertz  range,  transferring  enough 


radio  frequency  voltage  into  the  pre-amp  casing  to  influence  the 
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circuit.  The  rectified  voltage  effectively  overloaded  the  first  stage 
amplifier,  overriding  the  weak  return  signal  from  the  array  and 
destroying  the  data.  This  effect  went  undetected  in  previous  use 


This  thesis  also  modeled  the  planar  acoustic  array  as  an 
optical  aperture,  obtaining  the  Fourier  transform  of  the  far  field 
'diffraction'  pattern  via  computer  calculations.  The  technique  rests 
on  the  fact  that 

.  .  the  Fraunhofer  diffraction  pattern  is,  the  Fourier  transform 

of  the  field  distribution  across  the  aperture  .  .  [Ref.  18,  p.  412] 
The  quote  applies  to  electromagnetic  radiation,  but  Dr.  waiters 
correctly  construed  the  electromagnetic  amplitude  as  an  acoustic 
pressure,  obtaining  the  following  equation; 

(6)  P(kx,ky)  =  1/iR  Hgperture  P(x>y)  exP  (kxx  +  kyy)l  dxdy  i 

-  P(kx,ky)  is  the  pressure  distribution  in  the  image  plane, 

-  kx  and  ky  are  the  spatial  frequencies  along  the  x  and  y  axes  of 

the  aperture  respectively;  kj=  k*J/R,  where  J=  x  or  y;  k=27i/\ , 

-  P(x,y)  is  the  pressure  amplitude  at  the  aperture, 

-  exp  {i  (kxx  kyy)}  is  the  spatial  variation  of  the  aperture  and 

-  R  is  the  distance  to  the  image  plane;  i  accounts  for  temporal 

phase  distributions  across  the  aperture. 


The  model  can  be  thought  of  as  a  plane  coherent  acoustic  wave 


passing  through  the  aperture,  which  produces  a  diffraction  pattern 
(the  Fourier  transform)  in  the  far  field  (the  image  plane)  by  means 
of  the  Huygens-Fresnel  principle.  From  this  pattern  it  was  possible 
to  reasonably  predict  the  size  of  the  main  lobe  and  the  approximate 
angles  and  intensities  of  the  side  lobes.  This  modeling  was  done  for 
both  the  square  array  and  the  more  densely  packed  hexagonal  array. 
The  results  are  presented  in  this  chapter. 


A.  HARDWARE 


The  computer  used  in  this  phase  was  a  Hewlett  Packard  (HP) 
217  with  an  INFOTEK  Basic  compiler,  identical  to  the  computer  used 
in  the  echosounder  arrangement.  The  memory  requirements  for  this 
program,  however,  were  tremendous,  requiring  over  4  megabytes  of 
memory  when  using  the  512  x  512  grids.  This  lead  to  insertion  of  2 
additional  memory  boards  into  the  computer,  (one  from  INFOTEK  and 
one  from  EVENTIDE),  bringing  the  hard  memory  to  over  5  megabytes. 
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B.  SOFTWARE 


The  program  [Appendix  B]  utilizes  a  Fast  Fourier  Transform 
(FFT)  routine  supplied  by  INFOTEK.  The  idea  is  to  enter  the  aperture 
(array)  as  a  small  series  of  points  on  a  huge  grid.  Typically  the  grid 
size  would  be  256  x  256  points  square,  and  occassionally  512  x  512. 
The  array  takes  up  only  a  small  center  portion  of  this  grid, 
approximately  40  x  40  grid  points.  The  larger  the  grid  field  in 
relation  to  the  input  aperture,  the  better  the  resolution,  provided  the 
input  aperture  is  not  so  small  as  to  obscure  the  finer  details. 

The  aperture  area  was  entered  digitally  as  a  series  of  'ones’, 
signifying  a  uniform  coherent  plane  wave  at  the  aperture,  with  the 
remainder  of  the  field  entered  as  zeroes.  The  FFT  routine  then 
computed  the  far  field  Fourier  transform  of  the  array,  and  plotted 
the  output  intensity  in  relation  to  a  gray  scale.  The  gray  scale  was 
adjusted  to  encompass  3  decades  of  intensity. 


C.  RESULTS  AND  ANALYSIS 


A  sample  of  the  input  apertures,  both  square  and 
hexagonal,  are  shown  in  Figures  19  and  20.  Generally,  the  radius 
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of  each  element  is  only  4  grid  units,  with  a  background  field  of  256 
x  256  grid  units.  These  figures  have  been  enlarged  to  show  detail. 
The  output  Fast  Fourier  Transform  for  each  is  shown  in  Figures  21 
and  22  respectively.  The  center  portion  of  the  FFT  for  each  array  is 
the  lowest  spatial  frequency  for  the  array,  and  corresponds  to 
the  main  lobe.  Rings  further  out  correspond  to  the  side  lobes  at 
larger  and  larger  angles.  Rings  far  from  the  main  lobe  are  not 
physical  except  at  very  high  frequencies.  Interestingly,  as  the 
frequency  is  increased,  more  and  more  side  lobes  are  forced  into  the 
forward  hemisphere. 

The  program  calculates  the  beam  structure  of  a  single  piston 
driver  on  an  infinite  plane.  The  actual  source  is  not  a  piston  driver 
(but  the  approximation  is  close).  Huygens  wavelettes  propagating 
backwards  from  the  end  of  the  horns  reflect  from  the  support  plane, 
increasing  the  energy  of  the  actual  lobe  field.  Also,  the  physical 
side  lobes  can  extend  no  further  than  90°  from  the  axis  of  the  array 
(due  to  surface  absorption  and  reflection)  but  this  limitation  is 
ignored  by  the  program.  In  addition,  an  obliquity  factor  was  not 


included.  These  limitations  affect  the  graphical  output,  but  do  not 
significantly  affect  the  following  calculations. 

The  angles  of  the  output  beam  pattern  rings  depend  on  the 
frequency  of  interest  {5  kHz),  the  dimensions  of  the  actual  array  and 
the  input/output  dimension  of  the  aperture.  To  calculate 

information  from  these  printouts  it  is  necessary  to  use  optical 
aperture  equations  from  Hecht  and  Zajac  [Ref.  18].  The  half  angle 
formula  for  a  circular  aperture  is: 

i  (7)  sin  0 «  q/R  =  1 .22  U2a  where, 

I 

\  -  q  is  the  distance  from  axis  to  side  lobe  of  interest, 

|  -  a  is  the  radius  of  aperture, 

|  -  L  is  the  wavelength, 

i 

-  R  is  the  distance  to  the  far  field, 

-  0  is  the  half  angle  or  cone  angle  from  the  axis  to  q. 

I 

I 

The  hexagonal  array  is  not  quite  circular,  but  in  the  limit  it 
approaches  a  circle.  In  this  case  the  constant  1.22  is  simply 
replaced  by  F,  a  constant  to  be  determined  by  experiment. 

By  using  Brigham’s  text  [Ref.  19]  it  is  possible  to  convert 
equation  (7)  into  the  form: 

(8)  0  =  F  ARCSIN  (Ln/NAT)  where, 
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-  0  is  the  angular  measure  of  the  side  lobe  from  the  main  lobe, 

-  L  is  the  wavelength  of  interest  =  6.86  cm, 

-  NAT  is  the  total  width  of  the  sampled  grid  =  (ND/X), 

-  N  is  the  width  of  the  input  window  (grid  units;  1 28,256  or  512), 

-  D  is  the  actual  dimension  of  the  physical  array  (cm), 

-  X  is  the  dimension  of  the  input  array  in  grid  units, 

-  n  is  the  output  number  of  grid  units  to  the  side  lobe  of 

interest=(rz/w), 

-  r  is  the  distance  to  the  output  side  lobe  of  interest  (cm), 

{NOTE:  X  and  D  must  correspond;  r  and  n  must  correspond}, 

-  z  is  the  width  of  the  output  window  in  grid  units  (<  N/2), 

-  w  is  the  width  of  the  input  and  output  windows  in  centimeters. 

As  an  example,  consider  Figure  23.  For  this  example,  D  = 
38.735  cm  (measured  between  the  outer  edges  of  the  farthest  2 
speakers),  w  =  10.3  cm,  N  =  256  grid  points  and  z  =  128  grid  points. 
To  calculate  the  angular  extent  of  the  main  lobe  in  real  space,  first 
measure  the  distance  (in  cm)  from  the  corner  of  the  output  quadrant 
to  the  edge  of  the  main  lobe,  along  the  horizontal  axis  (Fig.  23). 
Here,  r  =  0.7  cm  and  X  =  40  grid  points.  Therefore,  equation  (8) 
becomes: 

(6. 86)(0.7)(1 28/10.3)  cm 

(9)  0  =  F  x  ARCSIN  .  =  F  x  13.9°  . 

(256)(38.735)/40  cm 
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The  half-angle  of  the  actual  array  is  only  12°  as  measured  on 
Figure  8.  This  means  that  F  *  0.86  or  1/F  »  1.16  which  is  very  close 
to  1.22.  Using  the  value  F  =  1.16,  then  the  90°  point  can  be 
calculated  for  the  output  field.  In  this  case, 

(256)(38.735) 

(10)  Sin  {(1.1 6)(90°)}  .  =  r  =  2.81  cm, 

(6.86)(128/10.3)(40) 
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which  means  everything  outside  the  corner  radius  of  2.81  cm  can  be 
ignoredas  it  is  non  physical. 

The  output  FFT's  were  correlated  to  the  actual  lobe  patterns  of 
Figures  7  and  8,  the  apex  and  side  mounted  patterns  for  the 
hexagonal  array,  and  the  results  are  plotted  in  Figures  24  and  25 
and  tabulated  in  Table  1.  The  amplitudes  were  normalized  by 
assigning  a  value  of  unity  to  the  main  lobe;  the  angles  are  measured 
counter-clockwise  from  the  standard  X  axis;  the  amplitude  scale  is 
in  decibel  units. 

The  measured  amplitude  was  calculated  using  INVLOG  (dB/20), 
where  *dB'  is  the  normalized  decibel  level  estimated  from  Figures  7 
and  8.  The  computed  amplitudes  relied  on  the  formula  20  LOG 
"Q/Max"  [see  Appendix  B,  lines  1390  and  2380  for  definitions  of  Q 
and  Max], 

The  angles  were  computed  using  a  simpler  form  of  equation  (8); 

(11)  0  =  ArcSin  {  (K-1)  LVNAt  },  where  K  is  an  index  from 

1  to  N/2,  and  essentially  allows  the  computation  of  the  angle  at 
each  grid  point  chosen.  The  measured  angles  were  taken  directly 


from  Figures  7  and  8. 
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Fig.  24.  Polar  Plot  of  Apex  Mounted  Hexagonal  Array 


. . . 


Fig.  25.  Polar  Plot  of  Side  Mounted  Hexagonal  Array 
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TABLE  1 .  COMPARISONS  BETWEEN  ACTUAL  LOBE  PATTERNS  AND 
COMPUTED  LOBE  PATTERNS. 


APEX  MOUNT  (  Fig.  7  vs.  24  ) 


Measured 


Computed 


0 

amDlitude 

0 

amDlil 

1st  min 

1  1° 

.025 

12° 

.002 

1 st  max 

1  6° 

.126 

17° 

.128 

2nd  min 

21  0 

.056 

23° 

.007 

2nd  max 

26° 

.095 

28° 

.108 

3rd  min 

33° 

>.01 

o 

03 

CO 

.006 

3rd  max 

o 

CM 

.056 

41° 

.027 

4th  min 

51° 

.018 

46° 

.007 

4th  max 

64° 

.056 

o 

CM 

CD 

.148 

Measured 


SIDE  MOUNT  (  Fig.  8  vs.  25  ) 
lured  Computed 


0 

amDlitude 

0 

amoiit 

1st  min 

12° 

.050 

14° 

.048 

1st  max 

18° 

.158 

19° 

.151 

2nd  min 

27° 

>.01 

28° 

.01 

2nd  max 

o 

V3 

CO 

.039 

36° 

.030 

3rd  min 

52° 

>.01 

o 

CD 

.024 

3rd  max 

61 0 

.025 

56° 

.029 

% 
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The  side  mounted  section  corresponded  to  measurements  along 
the  X  axis,  and  the  apex  to  measurements  along  the  Y  axis.  The  error 
of  the  measurements  was  grossly  estimated  as  approximately  10% 
for  the  minimums,  as  they  tended  to  be  obscured,  and  3%  or  less  for 
the  maximums.  Some  values  exceed  this  amount  considerably  and 
require  further  investigation. 

Generally,  however,  the  tabulated  values  agree  strongly  and 
suggest  that  this  technique  has  considerable  credibility  at 
predicting  side  lobe  structures  of  untested  arrays. 


This  thesis  demonstrated  that  increasing  acoustic  planar  array 
speaker  density  from  25  speakers/1452  cm2  (1.72  x  102 
speakers/cm2)  to  19  speakers/943  cm2  (2.02  x  10‘2  speakers/cm2) 
effectively  suppressed  the  lobe  spreading  from  a  70°  cone  angle  to 
less  than  a  45°  cone  angle,  with  a  reduction  in  side  lobe  intensity  of 
between  4  and  10  dB. 

Graphical  profiles  demonstrated  the  great  variability  and 
volatility  of  turbulence  structures  in  the  atmosphere  and  underline 
the  need  for  better  tools  to  measure  and  understand  this  structure. 
These  profiles  identify  gravity  waves,  thermal  plumes,  boundary 
layers  and  other  temperature  induced  structures.  Particulary 
important,  this  device  allowed  tentative  quantification  of  altitude, 
thickness  and  relative  intensity  of  mountain  boundary  layers  over  a 
potential  DOD  laser  site. 

Independent  verification  of  the  CT2  values  is  important  and  is 
currently  under  study.  It  will  further  enhance  the  utility  of  this 


device.  Due  to  the  strong  exponential  dependence  of  acoustic 


attenuation  (B)  with  increasing  frequency,  future  improvements 
should  also  include  lower  frequency  drivers  (2-3  kHz)  of  similarly 
high  efficiency,  to  improve  the  range  of  the  device. 

It  was  shown  that  acoustic  array  beam  patterns  could  be 
successfully  modeled  on  a  computer  by  means  of  a  Fast  Fourier 
Transform  algorithm  and  tested  for  optimal  design.  Future 

improvements  to  this  techniaue  include  a  method  to  stochastically 
simulate  turbulence  and  model  the  beam  structure  as  it  passes 
through  this  artificial  turbulence.  Calculation  of  main  lobe  energy 
versus  total  energy  for  various  arrays  should  also  be  studied. 
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APPENDIX  A 


List  of  Material  Sources 


Disclaimer:  These  sources  are  listed  only  in  case  exact 


duplication  of  this  thesis  is  required  at  some  future  date.  Other 


equally  competitive  sources  for  each  item  do  exist  and  their 


exclusion  is  not  a  reflection  of  product  quality  in  any  way. 


Barium  impregnated  foam: 
1 .  Barrier  Corp. 

Tigard,  Oregon 
(503)  639-4192 
POC:  Mark  Dove 


2.  Latchford  Corp. 

(714)  734-3660 
POC:  Tom  Pelligreeni 


Cylindrical  Enclosure: 

1.  Rubbermaid,  Inc. 

1147  Akron  Rd 
Wooster,  Ohio  44691 
(216)  264-6464 


2.  Traffic  Products  Co. 
Citrus  Heights,  Ca 
(916)  723-2958 
POC:  Michael  Fish 


Speakers: 

Motorola  Inc. 

4800  Alameda  Blvd  N.E. 
Albuquerque,  N.  M.  87113 
(505)  822-8801 


Temperature/humidity  Indicator 
Weathermeasure 
P.O.  Box  41039 
Sacramento,  CA  95841 
(916)  481-7565 
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APPENDIX  B 


Fast  Fourier  Transform  Program 


1  Program  plots  either  the  input  points,  the  FFT  or 
Poiar  Plots  of  a  19  element  hexagon  array  or  a  25 
1  element  square  array. 

1  This  program  is  interactive. 

1  Uiu  Dec.  19,  86/LRM  modified  in  OCT  87 
COM  fie, ! : 5 1 Z  )  ,Im(  I : 5 1 2  ) 

INTEGER  M  ,N 

DIM  Real (51 2 .512  )  ,  Imag(5l2 ,512 ) 
fi E PL  Amp (  103  >  ,Ang(  100) 

PRINT  ENTER  THE  GRID  SIZE  Ni  N=  512,256  OR  123" 
input  n 

PRINT 

i 

N  2  -  N  /  2 

M-L06< N >/L0G(  2 ) 


PRINT  "ENTER 
PRINT  "ENTER 


TO  SEE  THE  19  ELEMENT  HEXAGON. 
TO  SEE  THE  25  ELEMENT  SQUARE .  " 


TO  SEE  A  PLOT  OF  THE  INPUT  POINTS' 
TO  SEE  A  PLOT  OF  THE  FFT" 

TO  SEE  A  POLAR  PLOT  OF  INTENSITY" 


PRINT  "ENTER  ’2’  TO  SEE  THE  25  ELEMENT  SQUARE." 

INPUT  X 

PRINT 

PRINT 

PRINT  "ENTER  1  I’  TO  SEE  A  PLOT  OF  THE  INPUT  PO II 

PRINT  "ENTER  ’2’  TO  SEE  A  PLOT  OF  THE  FFT" 

PRINT  "ENTER  ’5’  TO  SEE  A  POLAR  PLOT  OF  INTENSI 

INPUT  A 

PRINT 

PRINT 

PRINT  "ENTER  THE  GRID  SIZE  OF  EACH  ELEMENT,  R . “ 
PRINT  "  USUALLY  R  =  4  OR  8. " 

INPUT  R 

PRINT 

PRINT 

PRINT  "ENTER  ’7’  TO  SEE  THE  FIRST  QUADRANT  ONLY 
FRINT  “  ELSE  ENTER  'O'." 

INPUT  C 

•  B=*  THE  SIZE  OF  THE  OUTPUT  WINDOW 
B»N2/2 


IF  X *  I  THEN 


F  =  2 . I »R*C0Si  .5235988  ) 

FOR  I =N2 -6 *R  TO  N2+6«R 
FOR  J-N2-5*R  TO  N2+6»R 

IF  <  <  I+4»R-N2  ) "2  +  (  J-N2  )  "2  >'R*R  THEN 
ReaK  I  , J  )=  I  .000 


FOR  K  = 1  TO  4 

Real( I+2»K*R  ,J  >=Real<  I  ,J  ) 
Real(  I+2*K*R-R  ,J+F  >=Real( I , J ) 
Real ( 1+2 ♦K*R-R  ,J-F  >=Real( I ,J ) 
NEXT  K 

FOR  K-1  TO  3 

Real(  I +2  *K  *R  ,J-2»F  >=Real( I ,J ) 
Real( I+2»K*R  ,  J+2*F  )  =  Real( I . J ) 

NEXT  K 


END  IF 
NEXT  J 
ME  a  I 

i 

END  IF 

i 

IF  X  =  2  THEN 

i 

K-N2+5»R 
L=NZ-5*R 
FOR  1=1  TO  N 
FOR  J-l  TO  N 

IF  I  < K  AND  I  >L  AND  J-'K  AND  J>L  THEN 
Cy=N2-4*R 
FOR  P=l  TO  5 
Cx  =N2 -4  »R 
FOR  0=1  TO  5 

IF  KCx  +  R  AND  I.Ca-R  AND  J'-Cy-R  AND  J,Cy+R  THEN 
IF  R*R  ( ( I-Cx  )"2  +  ( J-Cy  )"2  )  THEN 
Real( J  ,1  )=1 .000 
END  IF 
END  IF 
C*=Ca+2»R 
NEXT  0 
Cy=Cy  +  2  »R 
NEXT  P 
END  IF 
NEXT  J 
NEXT  I 

i 

ENO  IF 

PRINTER  IS  701 

i 

1  This  section  plots  the  input  points  only 
I 

IF  A= 1  THEN  GOTO  2240 

i 

•  THIS  SECTION  IS  THE  FAST  FOURIER  TRANSFORM 


T I = T I MEDATE 
FOR  1=1  TO  N 


9  T  0  Re(J)=Real<I,J) 

930  Im(  J  )=Imag( I  ,  J  ) 

990  NEXT  J 

1000  ' 

!  0  I  0  CALL  Fft(M) 

1 020  1 

1030  FOR  J=l  TO  N 

1040  Real( I  . J  )=Re( J  ) 

1090  Imag(I,J>=Im(J) 

1090  NEXT  J 

1070  NEXT  I 
1030  1 

1090  FGR  J=1  TO  N 
1 l 00  FGR  1=1  TO  N 

1110  Re<  I  >=Rea  1 (  I  .J  ) 

1100  I n ( I  )  = Imag( I  ,  J  ) 

1 1 30  NEXT  I 

1140  i 

1 150  CALL  Fft(M) 

115  0  ' 

1 1 70  FOR  1=1  TO  N 

I  1 80  Real ( I  , J  )=Re<  I  ) 

1  i  90  I  mag  (  I  ,  J  !=■!«(  I  ) 

1200  NEXT  I 

1210  NEXT  J 
1220  T  2  —  T IMEOATE 

1230  '  PRINT  INT( 10O(T2-T1  )/B0  )/ I  00  ,  "MINUTES ” 

1240  i 

1250  IF  A=2  THEN  GOTO  1 9S0 
I  2  G  0  1 

1270  1  This  section  graphs  the  polar  plot  of  the  FFT. 

12  30  i 

1230  B I =N2/2 
1300  82- 1 

1310  1  0=  The  actual  dimension  of  the  array(cm) 

1320  0=38.735 

I330  P m a >.  =  (  1  +  INT(  D*N/ (  68 . 6  *R  >  )  > 

1 340  ' 

1350  1  The  loop  does  the  side  and  ape.-,  mounts 

1 3E0  i 

1373  Rl=Real<  1,1) 
l 380  Ia=Imag( 1,1) 

1350  Na a  =  5QR ( R 1 *R 1  + 1  a » I  a > 

1400  1 

14  10  FOR  K  = 1  TO  2 

14  20  i 

1430  Ic  X = 1  AND  K“ 1  THEN 

1440  PRINT  “  HEXAGONAL,  ARRAY  ,  SIDE  MOUNTED 

1450  END  IF 

1460  IF  X  =  1  AND  K  =  2  THEN 

1470  PRINT  '  HEXAGONAL  ARRAY,  APEX  MOUNTED 


END  IF 

IF  X  =  2  AND  K  =  1  THEN 
PRINT  “ 

END  IF 


SQUARE  ARRAY  SIDE  COUNTED" 


1  This  section  conputes  the  ape*  square  array 

i 

IF  X  =2  AND  K  =  2  THEN 

PRINT  "  SQUARE  ARRAY,  APEX  MOUNTED" 

1  The  'L'loop  does  only  the  ape-  mounted  square 
FOR  L  = 1  TO  N2 
R  1  =Rea 1 ( L  ,L  ) 

G=Imag( L  ,L ) 

Q=SQR( R1 »R1+G*G ) 

IF  L  ■  P  m  a  a  THEN  GOTO  1 660 
Amp  (  L  )=Q/'Max 

Ang(L  i=ASN( ( (L-l  )  »63 . 6  * R  )/(D*N)  )  *  I  90/ P I 
i  PRINT  L  ,Q/Ma.<  ,ASN(  (  <  L—  1  )*68.6*R  )/(D*N  )  i  *  1  SO  /  3  .  1415927 
NEXT  L 

CALL  Polarplt(Amp(  »  )  ,Ang(  *  )  ,Pma.«.  ) 

IF  X  =  2  THEN  GOTO  2590 
END  IF 


FOR  1=1  TO  B 1 
FOR  J=l  TO  B2 
IF  B 1  =  t  THEN  P  =  J 
IF  82= 1  THEN  P=I 
R  1  =Rea 1 ( I  ,J  ) 

G= I mag( I  ,  J  ) 

G=SQR( R I »R 1 +G»G  ) 

IF  P,Pma*  THEN  GOTO  1840 
Amp( P  )=Q/Ma* 

Ang(P  )=ASN< <(P-1  )  *G8 . 6 • R  )  / ( 0 »N  )  )  • 1 80 /P I 
i  PRINT  P  ,Q/Max  ,ASN( < (P-1  )»S8.G*R  >/ ( 0*N  )  >*180/3. 1415827 
NEXT  J 
NEXT  I 

i 

CALL  Polarplt(Amp( ♦  )  ,  Ang<  *  )  ,Pma*  > 

i 

B2=B  1 
B 1  -1 
NEXT  K 

i 

IF  A=3  THEN  GOTO  2590 

i 

1  This  section  centers  the  F FT  onto  the  output  grid. 


IF  07  THEN 
GCLEAR 
G I N I T 


USE  ONLY  TO  GET  ONE  QUADRANT  OF  THE  GRID 


bnAF  H I CS  ON 
VIEWPORT  30  ,  I  30  ,0 . 103 
WINDOW  l,N2,!,N2 

6  I  =n:-i 
62=0 

GOTO  2340 
END  IF 

i 

1  *  r :  s  section  centers  the  output 


FOR  1  =  1  TO  N 2 -  1 
FOR  J=l  TO  N2 - 1 
R 1  =  R  e  3  I  (  I  ,  J  ) 

I  a  =  I  n  a  g  <  I  ,  J  ) 

Rea  1 <  N2  + 1 -  1  .NZ+J-1  )  =  R1 
R  e  a  1  <  N  2  - 1  ,  N  2  J  -  I  )  =  R 1 
Rea  1  <  TJ 2  —  I  ,N2-J  )=R1 
Rea  1 ( N2  + 1 -  I  ,N2-J  >  =  R1 
I  m  a  g ( N  2 -  I  ,N2-J  )  =  I  a 
Inag( N2  +  I-I  ,N2-J  >  = I  a 
I  nnag (  N2  ♦  I  -  I  , NZ+J-l  )=Ia 
I  ci  a  g  (  N  2  - 1  ,  N  2  +  J  - 1  >  ■=  I  a 
NEXT  J 
NEXT  I 


6CLEAR 

5INIT 

GRAPHICS  ON 

i 

VIEWPORT  30 , 1 30  ,0 , 1  00 
WINDOW  N2-B  ,N2  +  B  ,N2-6 ,N2  +  0 


FOR  I  =N2 -B I  TO  N2+62 
FOR  J  =  N 2 - B 1  TO  N2  +  B2 
R 1 =Rea 1 < I  ,J  ) 

Ia=Imag<  I  ,  J  ) 

Q  =  SQR(Rl*Ri  +  Ia»Ia  ) 

IF  A= 1  THEN 

Z  =  <  I2»Q+1  ) / 1  6  '  USE  ONLY  FOR  PLOTTING  INPUT  POINTS 

END  IF 
IF  A  =  2  THEN 

Z=  1  +LGT(  Q/Ma.*.  ) /G  1  GNLY  FOR  DECADE  INTENSITY  Sr'rEAD 
END  IF 

IF  Z-0  THEN  Z=0 
IF  Z' 1  THEN  Z= I 
AREA  INTENSITY  Z  ,2  ,Z 
MOVE  I  ,J 

RECTANGLE  1,1, FILL 


NEaT  I 

I 

PRINT  "The  sice  of  the  output  window  is"  , B  , " X "  ,8 
PRINT  "The  sice  of  the  input  window  is",N,"X",N 
PRINT  "The  radius  of  each  element  is  " ,R 

i 

1  Pause  statements  allow  insertion  of  parameters  w/o 
1  redoing  the  FFT  calculations. 

PRINTER  15  CRT 

PAUSE 

FmUSE 

END 

i 

1  This  section  is  the  FFT  subroutine  (from  INFOTEK). 

i 

5 UB  F  f  t  <  INTEGER  M) 

CON  Fe(  I  :5 I  0  )  ,  Im<  I  : 5  I  2 ) 

'  FFT  SUBROUTINE  FROM  INFOTEK  29  AUG ,  1985: 
i  USE  COMMON  FOR  FASTER  EXECUTION  3.09  sec  vs  3.59 
INTEGER  I,J,K,L,N,Ip,Le,Lel 
N=2  ‘  M 
J  =  1 

FOR  1=1  TO  N-l 
IF  I-J  THEN 
T  =  Re < J  ) 

Re ( J  >  =  Pe (  I  ) 

Re*  I  )  =  T 
T  =  I  m  (  J  ) 

I m(  J  )  =  I m(  I  ) 

I  m  (  I  >  =  T 
END  IF 
K  =N  DIU  2 
WHILE  K<J 
J  =  J-K 
K=K  DIU  2 
END  WHILE 
J  =  J  +K 
NEXT  I 
Le=  I 

FOR  L= I  TO  M 
Lei =Le 
Le=Le+Le 
Ure=  I 
Uim  =  0 
Ang  =  P I / L e  t 
Ur e  =  C0S ( Ang  > 

Uim=5IN(Ang) 

FOR  J  =  I  TO  Lei 

FOR  I=J  TO  N  STEP  Le 
I p  =  I +  Le  I 


,'.V7 


i  4®! 


I^e-fiei  Ip ) »  Ur  e - 1 n (  Ip )*Uin 
Tir  =  Re(  Ip  ) «U im  + I m(  Ip  )*Ure 
Ret  Ip  )  =  Ret  I  )-Tre 
I  r*  (  Ip  )  =  I  n  (  I  )  -  T  i  m 
Re  <  I  )  =  Re i  I  >  +  Tre 
I  r  f  I  )  =  I  n  (  I  )tT  in 
NEiT  I 

7-iJre*Ure-Ui'vi',Wir'i 
L  : "  -ijr e  *U  X  m  +  U  l  m  *  Ure 
Ure  =  T 
T  J 
r  l 

:  a  1  s  =  1  .  /  N 
Re  =  •;  Upscale  )*Re 
I m=  i Np sea  1 e  > *  I ^ 


1  '’'us  section  is  the  polar  plot  subroutine  written 

1  by  F  a u 1  Davison. 

I 

SuB  Fo i arp 1 1  <  Amp  i  *  i  ,Ang(  *  )  ,Pna>.  ) 

REM  THIS  MODULE  MAPS  THE  CALCULATED.  NORMALIZED 
REM  DSC  3  AFRAV  D  S  < •  >  ON  LOG  NORMAL  POLAR  PLOTS. 

I 

i  INITIALIZING  THE  PLOTTER  +++++++++++++»++++ 

i 

G  INI  T 

plotter  is  crt  internal" 

j F. APHlCi  uN 
6  C  L  E  A  R 

i 

UlEuPuRT  25,100,20,95 
ShOw1  -40,40,-40,40 

l 

_  I NE  TvPE  4 
R  i  n  g  =  1  0 

AxES  Ring.Ring.O.O.E.S.Z 

LINE  TYPE  3 

DEG 

LORG  7 
CSIZE  3 
FOR  R=0  TO  4 
LINE  T-FE  1 
MOLE  O  ,  R  »  R l ng 
l-6EL  P ‘Ring-40 
LINE  T  r'FE  3 
POP  Angle=0  TO  360 

PLOT  R*Rir.g»COS(  Angie  )  ,  R »  R  i  ng  »  S  I N  <  Angie  > 


£ 


3520  NEXT  Angle 

3530  PENUP 

3540  NEXT  R 

3550  FOR  Angle=30  TO  150  STEP  30 
3550  IF  Ang  1 e \ >  90  THEN 

3570  PLOT  40»COS< Angle  )  ,  40*S INI Ang 1 e >  ,-2 

3530  0RAU1  40»COS(  I  80+Angle  )  ,40«SIN(  I  80*-Angle  ) 

3530  END  IF 

3500  NEXT  Angle 

3610  PENUP 

3620  1 

3630  LINE  TYPE  1 
3640  OEb 

3650  FOR  P=Pna*  TO  1  STEP  -I 
3660  R=40+Z0*LGT< Amp( P  )  ) 

3670  IF  R>=0  THEN 

3G80  PLOT  R*COS(Ang(P  )  >,R»SIN(Ang(P)  ) 

3690  END  IF 

3700  NEXT  P 

3710  FOR  P=l  TO  Pmax 

3720  R=40+20»lGT< Anp( P >  ) 

3730  IF  R >  =0  THEN 

3740  PLOT  R*COS<Ang(P  )  )  ,-R»SIN(AngCP  )  ) 

3750  END  IF 

3760  NEXT  P 
3770  * 

3790  CLIP  OFF 
3790  1 

3800  LORG  2 
3310  MOVE  40,0 
3820  LABEL  "  X " 

3830  LORG  8 
3840  MOVE  -40,0 
3650  LABEL  ”-X  " 

3660  LORG  6 
3970  MOVE  0  ,-40 
3880  LABEL  "-Y" 

3890  1 

3900  CLIP  ON 
3910  1 

3920  DUMP  GRAPHICS  #701 
3930  GRAPHICS  OFF 
3940  PRINT  CHRS1 12 ) 

3950  SUBEND 
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